Constraints on Models for Proton-Proton Scattering from 

Multi-strange Baryon Data 
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Abstract 

The recent data on pp collisions at 158 GeV provide se- 
vere constraints on string models: These measurements 
allow for the first time to determine how color strings are 
formed in ultrarelativistic proton-proton collisions. 

1 Introduction 

Recently, the NA49 collaboration has published the 
rapidity spectra of p, A, S as well as the corresponding 
antibaryons in pp interactions at 158 GeV. These mea- 
surements provide new insight into the string formation 
process. In the string picture, high energy proton-proton 
collisions create "excitations" in form of strings, being 
one dimensional objects which decay into hadrons accord- 
ing to longitudinal phase space. This framework is well 
confirmed in low energy electron-positron annihilation j2] 
where the virtual photon decays into a quark-antiquark 
string which breaks up into mesons (M), baryons (-B) and 
antibaryons (B). An example of a q — q string fragment- 
ing into hadrons is shown in FigQ] Proton-proton col- 
lisions are more complicated due to the fact that even 
at 158 GeV proton-proton collisions are governed by soft 
physics, thus pQCD calculations can not be applied. And 
the mechanism of string formation is not clear, as will be 
discussed in the following. 
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Figure 1: e + e 
hadrons. 



qq. The q — q string fragments into 



• Longitudinal excitation (LE) models: UrQMD [3J, 
HIJING g], PYTHIA 0, FRITIOF (Ej; 

• Color exchange (CE) models: DPM 0, VENUS 0, 
QGS 0. 




One may distinguish two classes of string models: 



Figure 2: Two string formation mechanisms for pp colli- 
sions are presented: a) longitudinal excitation (LE). b) color 
exchange (CE). 



In the LE case the two colliding protons excite each other 
via a large transfer of momentum between projectile and 
target, FigEK- In contrast, the CE picture considers a 
color exchange between the incoming protons, leaving be- 
hind two octet states. Thus, a diquark from the projectile 
and a quark from the target, and vice versa, form color 
singlets. These are identified with strings, c.f. FigEh- 
The color exchange is a soft process. The transfer of mo- 
mentum is negligible. The final result, two quark-diquark 
strings with valence quarks being their ends, however, is 
quite similar. 

How are baryons and antibaryons produced? The eas- 
iest way to obtain baryons is to break the strings via 
quark-antiquark pair production close to the valence di- 
quark. Since the ingoing proton was composed of light 
quarks(qqq), the resulting baryon is of qqq or qqs type. 
Thus nucleons, As or Ss are formed. Since these baryons 
are produced at the string ends, they occur mainly 
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close to the projectile rapidity or target rapidity (leading 
baryons). 

Multi-strange baryons which consist of two or three 
strange quarks are produced near the quark end or the 
middle of the strings, via ss-ss production. Therefore the 
distributions of multi-strange baryons are peaked around 
central rapidity and the corresponding yields of multi- 
strange baryons and their antiparticles should be compa- 
rable. A closer look reveals an interesting phenomenon: 
Theoretically one finds the ratio of yields |9jj: 



0. 



Experimentally, however, S + s are less frequent than ex- 
pected. The ratio at midrapidity is 



1.2. 



0.44 ±0.08. 



The situation for fls is even more extreme: from string 
models one gets [2] 

n + /n- = i.6 - 1.9 

at midrapidity. From extrapolating A and 2 results (and 
from preliminary NA49 data) we expect 



0.5 - 0.8. 



This is a generic situation; It is impossible to get the 
fl ratio smaller than unity from those two types 

of string models. As addressed in this is due to the 
fact that the strings have a light quark (but not a strange 
quark) at the end, which disfavours multi-strange baryon 
production, and does not allow for fl production in the 
fragmentation region. 



a) 



b) 



Figure 3: Two chains of quark lines (b) representing a 
pair of strings (a) 
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Figure 4: Two pairs of strings (a) and the corresponding 
chains (b) 



where x is a function of energy and impact parameter, 
see fig. El Here, a dashed vertical line represents an 



prob( |)~ X prob( | i)~ -y 



2 Problems with the String Model 
Approach 

So is there something fundamentally wrong with string 
models? To answer this question let us consider some- 
what more in detail how string models are realized. One 
may present the particle production from strings via 
chains of quark lines [J] as shown in fig. It turns 
out that the two string picture is not enough to explain 
for example the large multiplicity fluctuations in proton- 
proton scattering at collider energies: more strings are 
needed, one adds therefore one or more pairs of quark- 
antiquark strings, as shown in fig. Q] The variables Xi 
refer to the longitudinal momentum fractions given to 
the string ends. Energy- momentum conservation implies 
X>i = 1- 

What is the probabilities for different string numbers? 
Here, Gribov-Regge theory comes at help, which tells us 
that the probability for a configuration with n elementary 
interactions is given as 

Prob(n interactions) = — - exp(— x), 



Figure 5: Probabilities for configurations with one and 
two elementary interactions (Pomerons), represented as 
dashed lines. 

elementary interaction (referred to as Pomeron). 

Now one identifies the elementary interactions 
(Pomerons) from Gribov-Regge theory with the pairs of 
strings (chains) in the string model, and one uses the 
above-mentioned probability for n Pomerons to be the 
probability for configurations with n string pairs. Unfor- 
tunately this is not at all consistent, for two reasons: 

1. Whereas in the string picture the first and the sub- 
sequent pairs are of different nature, in the Gribov- 
Regge model all the Pomerons are identical. 

2. Whereas in the string (chain) model the energy is 
properly shared among the strings, in the Grivov- 
Regge approach does not consider energy sharing at 
all (the x is a function of the total energy only) 

These problems have to be solved in order to make reli- 
able predictions. 
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3 Basic Ideas of Parton-Based 
Gribov-Regge Theory 

We consider a new approach called Parton-Based Gribov- 
Regge theory to solve the above-mentioned problems. 
Here we still use the language of Pomerons as in Gribov- 
Regge theory to calculate probabilities of collision con- 
figurations and the language of strings to treat particle 
production. Multiple interactions happen in parallel. An 
elementary interaction is referred to as Pomeron. The 
spectators of each proton form a remnant, see Fig. A 
Pomeron is finally identified with two strings, see Figd 
But we treat both aspects in a consistent fashion: In both 
cases energy sharing is considered in a rigorous way[10 , 
and in both cases all Pomerons are identical. This is 
the new feature of our approach, and the second point 
is exactly the solution to multistrange baryon problem 
mentioned above. In the following we discuss how to re- 
alize the collision configuration and particle production, 
respectively. 



der the condition that energy sharing is considered rigor- 
ously. Once a collision configuration is determined, not 
only the number of Pomerons but also the energy shar- 
ing among Pomerons is fixed. This is different from other 
string models. In the following we discuss how to realize 
it in a qualitative fashion. 

4.1 Reminder: some Elementary Quan- 
tum Mechanics 

Let us introduce some conventions. We denote elastic two 
body scattering amplitudes as T 2 ^ 2 and inelastic ampli- 
tudes corresponding to the production of some final state 
X as T 2 ^x (see figEl ) . As a direct consequence of uni- 
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Figure 6: Multiple elementary interactions (Pomerons) in 
NEXUS. 
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Figure 7: Multiple elementary interaction (Pomerons) in 
NEXUS (a). An elementary interaction is identified with a 
pair of strings. 



Figure 8: An elastic scattering amplitude T 2 _ 2 (left) and 
an inelastic amplitude T 2 _,x (right). 

tarity on has 2ImT 2 ^2 = Y, x ( T 2^x){T 2 ^x)* ■ The right 
hand side of this equation may be literally presented as 
a "cut diagram", where the diagram on one side of the 
cut is {T 2 ~,x) and on the other side (T 2 ^x)* , as shown 
in figEl • So the term "cut diagram" means nothing but 





X 



Figure 9: The expression J2x (^2->x)(22->x)* which may 
be represented as a "cut diagram". 

the square of an inelastic amplitude, summed over all fi- 
nal states, which is equal to twice the imaginary part of 
the elastic amplitude. Based on these considerations, we 
introduce simple graphical symbols, which will be very 
convenient when discussing multiple scattering, shown in 
fig. EH a vertical solid line represents an elastic am- 
plitude (multiplied by i, for convenience), and a vertical 
dashed line represents the mathematical expression re- 
lated to the above-mentioned cut diagram (divided by 
2s, for convenience). 



4 Collision Configuration 

Gribov-Regge theory is applied in NEXUS to calculate 
collision configurations. This calculation is performed un- 



4.2 Elementary Interactions 

Elementary nucleon-nucleon scattering can be considered 
as a straightforward generalization of photon-nucleon 
scattering: one has a hard parton-parton scattering in the 
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= iT,. 



~ 2s T 2^X T 2^X ~ 2s 2 Im7 2^2 



interactions to happen in parallel, where an interaction 
may be elastic or inelastic. The inelastic amplitude is the 
sum of all such contributions with at least one inelastic 
elementary interaction involved. To calculate cross sec- 



Figure 10: Conventions. 



middle, and parton evolutions in both directions toward 
the nucleons. We have a hard contribution Th ar d when 
the the first partons on both sides are valence quarks, 
a semi-hard contribution T som i when at least on one side 
there is a sea quark (being emitted from a soft Pomeron) , 
with a perturbative process happening in the middle of 
Pomeron, and finally we have a soft contribution when 
no perturbative process happens at all (see fig. ITTf . The 
total elementary elastic amplitude T2->2 is the sum of all 
contributions. 
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Figure 12: Inelastic scattering in pp. 

tions, we need to square the amplitude, which leads to 
many interference terms, as the one shown in fig. EI&), 
which represents interference between the first and the 
second diagram of fig. El Using the above notations, 
we may represent the left part of the diagram as a cut 
diagram, conveniently plotted as a dashed line, see fig. 
Elb) . The amplitude squared is now the sum over many 
such terms represented by solid and dashed lines. 




(b) 



Figure 11: The elastic amplitude T 2 ^2- 

The cut-off of virtuality between perturbative and non- 
perturbative processes, Qq, is independent of the beam 
energy. The perturbative part of Pomeron is strictly 
based on the standard DGLAP evolution with ordered 
parton virtualities in the ladder diagram and the non- 
perturbative part is based on Reggeon parameterization. 
The semihard Pomeron is a convolution of hard and soft. 
The total elementary elastic amplitude Ti^2 is the sum 
of all these terms. Thus we have a smooth transition from 
soft to hard physics: at low energies the soft contribution 
dominates, at high energies the hard and semi-hard ones, 
at intermediate energies (that is where experiments are 
performed presently) all contributions are important. 

The multiple scattering theory will be based on these 
elementary interactions, the corresponding elastic ampli- 
tude T2_>2 and the corresponding cut diagram, both be- 
ing represented graphically by a solid and a dashed ver- 
tical line. We also refer to the solid line as Pomeron, to 
the dashed line as cut Pomeron. 

4.3 Multiple Scattering 

We first consider inelastic proton-proton scattering, see 
fig. El We imagine an arbitrary number of elementary 



Figure 13: Inelastic scattering in pp. a) An interference 
term of cross section, b) Represented with our simplified 
notations. 

When squaring an amplitude being a sum of many 
terms, not all of the terms interfere - only those which 
correspond to the same final state. For example, a sin- 
gle inelastic interaction does not interfere with a double 
inelastic interaction, whereas all the contributions with 
exactly one inelastic interaction interfere. So consider- 
ing a squared amplitude, one may group terms together 
representing the same final state. In our pictorial lan- 
guage, this means that all diagrams with one dashed line, 
representing the same final state, may be considered to 
form a class, characterized by m = 1 - one dashed line 
( one cut Pomeron) - and the light cone momenta x + 
and x~ attached to the dashed line (defining energy and 
momentum of the Pomeron). In fig. El we show sev- 
eral diagrams belonging to this class, in fig. El we show 
the diagrams belonging to the class of two inelastic in- 
teractions, characterized by m = 2 and four light-cone 
momenta x±, x^ , x\ , x 2 . Generalizing these consider- 
ations, we may group all contributions with m inelastic 
interactions (m dashed lines = m cut Pomerons) into a 
class characterized by the variable 

K {fn, Xi , Xi , • • • , x m , x m }. 
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Figure 14: Class of terms corresponding to one inelastic 
interaction. 




Figure 15: Class of terms corresponding to two inelastic 
interactions. 

We then sum all the terms in a class K, 

fl(K) = { all terms in class K } . 

The inelastic cross section is then simply a sum over 
classes, 




Q, depends implicitly on the energy squared s and the 
impact parameter b. The individual terms / d 2 bfl(K), 
represent partial cross sections, since they represent dis- 
tinct final states. They are referred to as topological cross 
sections. One can prove 

K 

which is a very important result justifying our interpre- 
tation of Q(K) to be a probability distribution for the 
configurations K . This provides also the basis for apply- 
ing Monte Carlo techniques. 

The function Q, is the basis of all applications of this 
formalism. It is the foundation not only for calculating 
(topological) cross sections, but also for particle produc- 
tion, thus providing a consistent formalism for all aspects 
of a nuclear collision. 

5 Particle Production 

Here we discuss how to produce hadrons with any given 
collision configuration determined according to the above 
section. As we discussed above, each cut Pomeron is iden- 
tified as a pair of strings. Similar to other string models, 
a Lorentz invariant string fragmentation procedure pro- 
vides a transformation from strings to final-state hadrons 
in NEXUS. However the string formation mechanism in 
NEXUS is very different from others because Pomerons 
are treated identically for calculating particle production. 



5.1 Hard Pomerons 

As discussed above, each elementary interaction has 
hard, soft and semihard contribution. So each Pomeron 
has certain probabilities to be of type soft, hard and 
semihard. Once a collision configuration is determined 
according to the above section, the number of each type 
of Pomerons is fixed. To give a proper description of 
deep inelastic scattering data, hard and some of semi- 
hard Pomerons are connected to the valance quarks of the 
hadron. In order to conserve the initial hadron baryon 
content and to keep the simple factorized structure with 
the leading logarithmic approximation of quantum Chro- 
modynamics, we associate a "quasi-spectator" antiquark 
(of the same flavour) to each valence quark interaction. 
So a hard Pomeron is a two-layer ladder diagram where 
one of them is hard with ordered virtualities and the 
other is soft. The hard ladder diagram gives a kinky 
string (each string segment has a constant velocity and 
finally is identified with those emitted perturbative par- 
tons) . The fraction of hard Pomerons in total Pomerons 
is very subdominant in average proton-proton collisions. 
At SPS energies hard and semihard Pomerons do not con- 
tribute at all, and even at collider energies the Pomerons 
connected to valence quarks are rare. Therefore, even 
though all types of Pomerons are included in the Monte 
Carlo, their contribution can be ignored in the following 
discussion. 

5.2 Soft Pomerons 

How to form strings from soft Pomerons? No matter 
whether single-Pomeron exchange or multiple-Pomeron 
exchange happens in a proton-proton scattering, all 
Pomerons are treated identically. Because of this, it is a 
natural idea to take quarks and antiquarks from the sea 
as string ends, because an arbitrary number of Pomerons 
may be involved, and valence quarks are not always avail- 
able to be string ends due to their limited amount. This 
point is different from the above-mentioned string mod- 
els, where all the string ends are valence quarks. Letting 
all the valence quarks stay in remnants, thus, string ends 
from cut Pomerons have complete flavour symmetry and 
produce particles and antiparticles in equal amounts. In 
the old version of NEXUS used in |Q] , a sea quark and an 
antiquark of the same flavour have been taken from the 
sea as string ends to keep flavour conservation. In the 
new version, NEXUS 3.0, the flavours of the string ends 
are independent. In order to compensate the flavour, 
whenever a quark or an antiquark is taken from the sea 
as a string end, a corresponding antiparticle from the sea 
is put in the remnant nearby. 

5.3 Remnant 

Remnants are a new object, compared to other string 
models. The partonic content of a remnant is very clear: 
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three valence quarks and the compensated partons. From 
this partonic structure one can estimate that its mass is of 
the order of proton mass. The 4- momentum of a remnant 
is that of an initial proton minus those being taken away 
by participant partons. So the 4-momenta of projectile 
remnant and target remnant are fixed once the collision 
configuration is given. The remnant mass can be calcu- 
lated according to its 4-momentum and the on-shell con- 
dition, however this value of mass is not reliable. During 
the configuration calculation one distributes energy of or- 
der 100 GeV to Pomerons and remnants where a proton 
mass is negligible. We justify the remnant mass with a 
distribution P(m 2 ) oc (m 2 )~ a , m 2 S (rn^^, x + s), where 
s is the squared energy at center mass system, m m ; n is 
the minimum mass of hadrons to be made from the rem- 
nant's quarks and antiquarks, and x + is the light-cone 
momentum fraction of the remnant which is determined 
in the collision configuration. Through fitting the data 
we determine the parameter a — 1.5. Remnants decay 
into hadrons according to n-body phase space[llj. 

5.4 Leading Order Discussion 

The configuration of leading order has only one cut 
Pomeron. The most simple and most frequent colli- 
sion configuration has two remnants and only one cut 
Pomeron represented by two q — q strings as in FigTlBti . 
Besides the three valence quarks, each remnant has ad- 
ditionally a quarks and an antiquark to compensate the 
flavour. 
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Figure 16: a) The most simple and frequent collision con- 
figuration has two remnants and only one cut Pomeron rep- 
resented by two q — q strings, b) One of the q string-ends can 
be replaced by a qq string-end. c) With the same probability, 
one of the q string-ends can be replaced by a qq string-end. 

The leg of a cut Pomeron may be of qqq type with 
small probability P qq , which means the corresponding 
string ends are a diquark and a quark. In this way we 
get quark-diquark (q-qq) strings from cut Pomerons. The 
qqq Pomeron leg has to be compensated by the three 
corresponding antiquarks in the remnant, as in FiglTfih. 
The (3q3q) remnant may decay into three mesons (3M) or 
a baryon and an anti-baryon (B+B). Since the 3M mode 
is favored by phase space, we neglect B+B production 
here. 



For symmetry reasons, the leg of a cut Pomeron is of 
qqq type with the same probability P qq . This yields a 
q — qq string and a (6q) remnant, as shown in Figllfib. 
The (6q) remnant decays into two baryons. Since q-qq 
strings and q — qq strings have the same probability to 
appear from cut Pomerons, baryons and antibaryons are 
produced equally. However, from remnant decay, baryon 
production is favored due to the initial valence quarks. 

6 Results 

Here, we will concentrate on baryon-antibaryon produc- 
tion, because there we obtain strikingly different re- 
sults compared to other models. However, we carefully 
checked as well mesons - essentially pion and kaon rapid- 
ity and transverse momentum spectra, where the results 
are quite close to what we obtained earlier with NEXUS 
2 or VENUS. 

FigEI depicts the rapidity spectra of baryons and an- 
tibaryons from NEXUS 3.0 with P qq = 0.02 (solid lines). 
As a comparison, we also show the preliminary data from 
the NA49 experiment (points). The contributions of 
projectile remnants (dashed lines), target remnants (dot- 
ted lines) and cut Pomerons (dashed dotted lines) to par- 
ticle production are also shown respectively in FigEI 
FigEI demonstrates that NEXUS 3.0 describes reason- 
ably the rapidity spectra of baryons and antibaryons in 
pp collision at 158 GeV. 

We also provide the particle yields at midrapidity, y G 
(2/cm - 0.5 , y cm + 0.5), from NEXUS 3.0, Pythia 6.2 and 
compare them to data in table 



yield 


NEXUS 3.0 


Pythia 6.2 


NA49 data 


P 


9.12x10^ 


4.85x10^ 


9.28x10^ 


P 


2.00x10^ 


1.64x10^ 


2.05x10^ 


A 


1.61X10" 2 


7.53xl0- 3 


1.79xl0~ 2 


A 


5.85xl0~ 3 


4.02xl0~ 3 


5.57xl0^ 3 




8.08xl0~ 4 


2.53xl0~ 4 


7.08xl0~ 4 




4.71 xl0~ 4 


2.20xl0~ 4 


3.12xl0~ 4 


fi- 


2.79xl0~ 5 


2.33xl0- y 






2.16xl0~ 5 


2.94xl0~ 6 





Table 1: Particle yields at midrapidity in pp collisions at 
158 GeV. 

From NEXUS 3.0, we get the ratios at midrapidity 

= o.58, n + /rr = o.77. 

In conclusion, it seems that old string models fail 
to reproduce the experimental 5 /S~ and anticipated 
f2 + /£l~ ratio so far. The string formation mechanism as 
employed in NEXUS 3.0 is able to reproduce the exper- 
imental data nicely. The rapidity distributions of multi- 
strange baryons as well as As and protons can be under- 
stood. The main point is the fact that the final result of 
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a proton-proton scattering is a system of projectile and 
target remnant and in addition (at least) one Pomeron 
represented by two strings. At SPS energy, the soft 
Pomerons dominate. In general the soft Pomerons are 
vacuum excitation and produce particles and antiparti- 
cles equally because their string ends are sea quarks. The 
valence quarks stay in remnants and favour the baryon 
production as compared to antibaryon production. 
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Figure 17: Rapidity spectra of baryons and antibaryons 
calculated from NEXUS 3.0 (projectile remnant contribu- 
tion: dashed lines; target remnant contribution: dotted 
lines;Pomeron contribution: dashed dotted lines; sum: solid 
lines) and NA49 experiment (points). 
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